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Abstract 

The physical consistency of the match of piecewise-C'' metrics is 
discussed. The mathematical theory of gravitational discontinuity hy- 
persurfaces is generalized to cover the match of regularly discontinuous 
metrics. The mean- value differential geometry framework on a hy- 
persurface is introduced, and corresponding compatibility conditions 
are deduced. Examples of generalized boundary layers, gravitational 
shock waves and thin shells are studied. 
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1 Introduction 

Is it possible to define weak solutions of the Einstein equations of class 
piecewise-C°, i.e. to generalize the compatibility conditions which replace 
the field equations on a singular hypersurface to the case when the metric is 
regularly discontinuous? 

To reach this goal would probably mean to define the most general class 
of regularly discontinuous weak solutions of the Einstein equations. It seems 
that this problem was never studied before in the literature. But, before we 
proceed, we need to discuss whether we are talking of something mathemat- 
ically and physically consistent or not. 

A fundamental concept of Riemannian geometry is that at any point of 
a submanifold there are coordinate choices for which the metric reduces to 
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the Minkowski flat metric. Clearly, if this choice is made on both sides of 
the discontinuity surface, any "jump" in the metric disappears. Thus, the 
metric discontinuity appears as a coordinate dependent concept, which is 
neither geometrically (nor physically) acceptable in the context of General 
Relativity. 

But we also have to consider that regularity of the global coordinates plays 
an important role in our approach, which is that of [T] and of the literature 
cited therein. In particular, since here the spacetime is only C^, we are led to 
considering {C^, piecewise C^) coordinate transformations. If the metric is 
discontinuous in some globally C° chart, it is in general impossible to obtain 
the vanishing of the metric jump on both sides of a hypersurface with a C° 
coordinate transformation (see section 2). Moeover in the following we are 
led in a natural way to considering coordinate transformations; the metric 
discontinuity is a tensor with respect to such coordinate changes! 

In other words the jump of the metric has a precise mathematical mean- 
ing, if we consider it in connection with global regular coordinates. 

In a well consolidated procedure, the assumption of continuity for the 
metric across a gravitational interface is usually taken for granted; however 
it follows from the limiting process of the thin sandwich modelization, in 
consequence of the hypothesis that the external derivatives of the metric 
are bounded [2]. Yet in this paper we are going to see that, even removing 
the assumption of continuity, it is still possible to define a generalized inner 
geometry of the discontinuity hypersurface; one thus can consistently find 
a corresponding generalized set of compatibility conditions, which obviously 
reduces to the usual ones when the continuity hypothesis is restored. 

Yet, which are the physical motivations to move to such generalization? 
Actually gravitational shock waves and thin shells are usually defined by the 
presence of singular curvature with a "delta" component concentrated on a 
hypersurface, situation which is well cast within the classic C° piecewise-C^ 
match of metrics [H [3]. 

We were originally led to consider solutions of class piecewise-C°, as pos- 
sible generalizations of shock waves and thin shells, by the sake of mathe- 
matical completeness, with the idea that phisical interpretation would follow. 
We actually found a reacher framework than the usual one, with some in- 
teresting new features (and even some rather undesiderable ones), which we 
display in this paper. 

There are two main theories in the literature for solutions of class C° 
piecewise-C^, i.e. that in terms of the second fundamental form (heuristic 
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theory, see e.g. [HE]) and that in terms of the curvature tensor-distribution 
(axiomatic theory, see e.g. [SI E]); such are equivalent through appropiate 
extensions (for a self-contained overview see e.g. [1]). 

The axiomatic theory appears to be inappropriate to the study of general- 
ized solutions, since the theory of distributions is basically linear. Even if we 
could in principle replace the discontinuous metric with its associated distri- 
bution , then it would be impossible to define, for example, replacements 
for the Christoffel symbols, since this would involve product of distributions, 
which, as it is generally believed, is impossible to define. In fact it was proved 
by Schwarz [^ that, under reasonable hypothesis, there can be no definition 
of commutative and associative operation on distributions which reduce to 
ordinary multiplication on integrable distributions (say on regular functions); 
thus in a word it is impossible to define product of distributions. 

Or is it? Colombeau [H O [10] developed a theory which apparently con- 
tradicts Schwarz's result. He introduced a very broad space of generalized 
functions, which extends the usual space of distributions, a subspace of which 
corresponds, in a certain sense (the correspondance is not 1 to 1), to usual 
distributions. Colombeau's formalism permits multiplication of generalized 
functions; but the contradiction with the impossibility theorem is only ap- 
parent, in fact Schwarz's hypothesis are violated, since the operation does 
not coincide with ordinary multiplication on regular functions nor with mul- 
tiplication of a regular function times a distribution (although it does at least 
for C°° functions). 

Such theory, however, does not fit in a natural way in general relativity, 
since it is impossible to define covariantly invariant geometrical objects; in 
fact Colombeau's space is not invariant for smooth coordinate transforma- 
tions, unless they are linear. Such difficulty, however, seems to have been 
overcome in subsequent adjustments of the theory, with the introduction of 
a richer mathematical framework [TTl [12] . so that the generalized functions 
current apparatus can be used in general relativity, and indeed it has been 
applied at least to the calculation of singular curvatures of the spacetimes 
of Kerr [13] , Reissner-Nordstrom [13] , and so-called cosmic-string spacetime 
[15] . In such literature Colombeau's theory is adapted to the handling of 
curvature when the metric has a singularity in the sense of functions, i.e. the 
ordinary curvature would explode, at a singular event-point or at a singular 
worldline. There seems to be no particular reason to forbid Colombeau's 
method also for defining the match of piecewise-C° regularly discontinuous 
metrics at a singular hypersurface; however, as far as the author is aware, no 
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attempt has been made yet to use it in this framework. 

The direct method we will introduce in the following sections, however, 
is so conceptually simple that we prefer not to experiment with Colombeau's 
generalized functions, which would instead mean introducing a far more com- 
plicated and unfamiliar mathematical apparatus. 

In this paper in fact we propose a new, generalized theory for regularly 
discontinuous solutions, covering also the match of piecewise-C° metrics. Our 
theory is heuristic, as it is constructed in a way similar to the heuristic the- 
ory of C°, piecewise-C^ solutions originated from the studies of Israel, but 
we completely avoid the traditional or projectional Gauss-Codazzi framework 
(which either does not include the lightlike case [HE], or needs a special adap- 
tation for it [TBI [I]) cind introduce what we called "mean- value differential 
geometry" framework, instead (see section 3). This is conceptually very sim- 
ple, and permits to construct in a natural way a generalized theory, where 
the main role (which used to be that of the jump of the secund fundamental 
form) is here played by the jump of the Christoffel symbols. The theory is an 
extension of the theory of gravitational discontinuity hypersurfaces we have 
studied in [I], to which it reduces when the metric is C°. Even if we should 
restrict to C° solution, by adding the traditional assumption of continuity for 
the metric, our theory would undoubtedly have at least the good qualities of 
not needing the timelike and the lightlike case to be distinguished (different 
from usual heuristic theory), and of just requiring C° continuity for the co- 
ordinates (different from the axiomatic theory). Moreover, it is completely 
cast in the framework of general coordinates of the ambient (glued) space- 
time, with no use of parametric equations of the hypersurface, nor of inner 
coordinates and holonomic 3-basis, which could be considered a good quality 
in some applications as well. 

Piecewise-C° weak solutions of the Einstein equations, as far as the au- 
thor is aware, have never been considered previously in the literature. They 
generalize the corresponding C° solutions, as examples in this paper show; 
however there is more. Apparently in fact the theory allows the propagation 
of free gravitational discontinuity at lower speed than the speed of light (sec- 
tion 8); or rather, we still have no general proof that the absence of stress 
energy concentrated on S should, in the timelike case, necessarily imply the 
degeneracy of a generalized solution to a boundary layer, although it does 
at least for a wide class of spherical matchs (see section 6). Moreover, non- 
simmetric stress-energy is allowed on the hypersurface (section 9), like e.g. in 
Einstein-Cartan dynamics. This possible link to classical unification theories 
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is surprising, since in our framework we have nothing similar to Einstein- 
Cartan torsion. We therefore see a lot of space for future investigation. 

2 Discontinuous metrics 

Let us suppose V4 an oriented differentiable manifold of dimension 4, of class 
(C°, piecewise C^), provided with a strictly hyperbolic metric of signature 
— and class piecewise-C°. Let C V4 be an open connected subset with 
compact closure. Let units be chosen in order to have the speed of light in 
empty space c = 1. Greek indices run from to 3. 

Let S C f2 be a regular hypersuperface of equation f{x) = 0; let f2+ and 

denote the subdomains distinguished by the sign of /. We suppose the 
metric and its first and second partial derivatives to be regularly discontin- 
uous on S in all charts of class C^n). Let / G n C^in\J:), and let 
second and third derivatives of / be regularly discontinuous on S. Finally, 
let ia = daf denote the gradient of /. 

Let the metric be a solution of the ordinary Einstein equations on each of 
the two domains and . In this situation S is the interface between two 
general relativistic spacetimes and it is called a (generalized) gravitational 
discontinuity hypersurface. 

In the following we will develope a theory to justify the introduction of 
suitable generalized compatibility conditions to replace the Einstein equa- 
tions on S (section 5); if such conditions are satisfied the match across the 
generalized gravitational hypersurface S will be called a generalized regularly 
discontinuous solution of the Einstein equations. 

Now let us briefly recall some basics notions on regularly discontinuous 
fields, which we will use as tools. In any case, for notation and terminology 
we refer to [1]. 

A field is said to be regularly discontinuous on S if its restrictions to 
the two subdomains $7+ and Q~ both have a finite limit for / — > 0; we 
denote such limits by ip~^ and ip~ , respectively. 

In this case the jump [ip] across S and its arithmetic mean value Tp are 
well defined on the hypersurface: 

W\ = 

Tp = (l/2)(^+ + v.-) 
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If is continuous across S, we obviously have: [ip] = 0, ip = ip. We also have 
the converse formulae: 




(2) 



(3) 



If a field ip is regularly discontinuous on S, its jump [ip] is sometimes called 
its discontinuity of order 0. 

The jump of a regularly discontinuous function has support on S, but in 
general, the partial derivative of the jump is well defined as the jump of the 
derivative of the function (see pTl [T8]). In particular, the derivative of the 
jump of a continuous field is not null, unless the field is also C^. 

Similarly, we define the partial derivative of the mean value as the mean 
value of the partial derivative. We can also use regular prolongations to 
extend, in a sense, the definition of Tp and [ip] to the whole domain f2. Thus 
they can be regarded as regular and derivable fields in fl, but their values 
(and those of their derivatives) are well defined only on E, while in Q\E 
they depend on the choice of the prolongation. For details on the method of 
regular prolongations see e.g. [T7l [TH]. 

We moreover define the covariant derivative of a field with support on S 
by means of the mean value Tpp'^ of the Christoffel symbols. For the jump of a 
regularly discontinuous vector, for example, with this definition one has that 
the jump of the covariant derivative is different than the covariant derivative 
of the jump. Thus, by definition, we have: 



and similarly for the jump of any regularly discontinuous tensor. 

Since the spacetime is only C°, we are led to considering (C°, piecewise 
C^) coordinate transformations, with regularly discontinuous first deriva- 
tives; the metric discontinuity [ga/s] is not a tensor with respect to such 
changes of coordinates. In fact we have: 



(4) 



and in consequence of ([3]): 



(5) 
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[gap] = [g, 



dx" dxf^ 



dx"' 
dx^ 



dx" 



where: 



i[g^ 



a'f3' 



dx^' 



+ ga'13' 



dxP' 

d^ 



(6) 



(7) 



We therefore can simulate all (C°, piecewise C^) coordinate changes by com- 
bining changes with metric gauge changes: 



[gap] ^ — ^ [gap] + qaP' 



dx^' 



+ Qa'P 



dx" 
dx" 



(8) 



which generalize usual gravitational gauge changes of the theory of (C°, piece- 
wise C^) solutions [Ij. 

Is it always possible to make [gap] vanish with an appropriate C° trans- 
formation? Clearly the answer is negative. In fact it suffices to consider the 
case when [gap] and (jap are both definite positive in a given chart to see that 
the equation obtained from by replacing the first hand side with has no 
solution for [dx" /dx"] and dx"' /dx". Thus the set of effective generalized 
gravitational discontinuity hyper surf aces is non empty. 

Moreover it will occur in many applications to have ^a ^ Therefore 
it will be often desiderable to work in the framework of (C^, piecewise C^) 
coordinate transformations, which preserve such condition. The metric dis- 
continuity is a tensor with respect to such changes of coordinates, but the 
jump of the Christoffel symbols, which appear to play a main role in the 
following, is not; we have in fact: 



[r 



aP 



[r. 



a'P' 



dx"' dx^' dx'^ 
^ dx" dxf^ dx""' 



+ 



d^x"' 



dx"dx^ 



dx" 
dx"' 



If the coordinates are and so is the form £„ we can write: 



-a'^P 



d^x"' 



(9) 



(10) 



_dx"dxf^ . 

where d"^ denotes the weak discontinuity of order 2 (see e.g. [I71[T8]). Thus on 
E we can generate all (C^, piecewise C^) transformations for [F] combining 
transformations (with respect to which F is a tensor) and Christoffel 
gauges transformations, i.e. of the kind: 
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[r„/]^[r./3l+CW (11) 

with some analogy with the case of C° metrics (where the main role is played 
by the first order metric discontinuity dg, see [1] section 3). 

In any case neither the mean value of the metric ^ or its jump [g\ now 
have null covariant derivatives. Consider in fact the identity V aQpp = in 
the domain fi"*"; from the limit / — > 0^, on S we have: 

d^gl - (r,/)+(7+ - {V^.'^ygl = (12) 

Here, with obvious meaning of the symbols, we denote: g'^p = {gpp)^ , g^p 
'g^, etc. Consequently on S, from (l2])i we have: 



9a9i3p + (l/2)9„ [c//3p[ - Ta/s^g^p - Tap''g^p+_ 

-(i/2)([r„/]^,^ + r„/[(7,,] + [r„/]^,^ + T^p^M+ (13) 
-(i/4)([r„/][^?,,] + [r,/]W) = o 

Similarly, from the limit / — > 0^ and from we also have on E: 

da9f3p - (1/2)9„[^/3pL- ^'^f^'^dup - T'ap''g„f3+_ 

+ (l/2)([r„/]^,^ + T^p^lgp,] + [r„/]^,^ + T^p'^M+ (14) 
-(l/4)([V][(7.p] + M[(7^,])=0 

From the sum of expressions ( |T3l) and ( |T^ we thus have: 

V^g^p = (l/4)([r,/][(7,,] + [r„/][(7^,]) (15) 
and, from difference: 

da[gi3p] = [Tafip] + [Tapf3] (16) 

From (1161) . ([3]), and from the definition of covariant derivative over S, we 
then have: 

^a[gpp] = [^al3'']g^p + [^ap'']gi3u (17) 

As for the jump and the mean value of the Christoffel symbols we have, from 
0: _ 

= {l/2){g'"'{dagfs^ + dpg^^-d^g^f^)+ . . 

+ il/A)[g--]idM + d^[g.^] - 9.M)} ^'^^ 
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and 

[r,/] = {i/2){r''mgpA + dp\g.a\-d,\g^p\)+ 

+ [g''"] (dagfsa + dpg^a - d^g^p) 

or, from ( |T5i) and ( |T7I) : 

3 Mean- value geometry on a hypersurface 

Let us consider a 4- vector regularly discontinuous on S, the jump and the 
mean value of which will work as a prototype of vectors with S as support. 
We have, by definition: 



[VpV^V''] = V4V„n - [r^„lV.y- + [V^z/-^] V„l^- (21) 
where [V«1/'"] = V aiV] + [Tau^W and where, again by definition, we have: 



v.^'^ = \{dAy^r + (r+).A'^(i^+)' + d,{V'Y + (r-).,'^(\/-)^} (22) 

Thus, from ([2]) we have: 



v,v- = v.r- + (i/4)[r,Al[n, (23) 

which, incidentally, is the same result we could get from the formal applica- 
tion of ([3]), wich actually can be apphed to covariant derivatives, provided 
one interpretes V = V. We therefore have: 

-[r,„iv.r - {i/W,A[^^xW']+ (24) 
+[r/3.iv„i^' + (i/4)[r^,-][r,Al[^^] 

and thus, by antisymmetrization: 

[V[^v„]n = Vi^v.jiv'^] + V[^[r„]..]F'^ + ^[r.[/][r„]/][v^] (25) 

Now, from the Ricci identity we have: [V[/3Va]V^°^] = [Ra(ip"VP] and then, by 
(ED: 
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[^[P^a]^] = [Rapp^V + Rcpp^iVP], (26) 

and thus from a well known identity which follows from ([3]) as a consequence 
our definition ([S]) for the covariant derivative on S, i.e. (see [I]): 

K/3p1 = V^[r,/]-V,M (27) 

we have that the commutator of the covariant derivatives of the jump of a 
generic regularly discontinuous vector obeys the following Ricci-like formula: 

V^p^^iV^] = {(1/2)7?.^/ - (l/4)[r.[/][r.]/]}[\/'']. (28) 

Not surprisingly, working in a similar way starting from V^VaV'^ and anti- 
symmetrizing, we find again: 

Vi^^jF"^ = {ii/2)R^pp^ - (i/4)[r,[/][r„]/]}F^ (29) 

in fact any given field with support on S can be considered, by the help of 
suitable prolongations, as the jump (or as the mean value of) some regularly 
discontinuous field. Thus, for any vector V with support on E, we can 
introduce the following mean-value geometry Ricci-like formula on E: 

{V[fsVa]W = (Rj^Up^VP; (30) 

where we have introduced the mean-value geometry curvature (-Re), defined 
by the following mean-value geometry first Gauss- Codazzi identity: 

{R^)app'' = Ra(3p" - (i/4)([r^,'^][r,/] - [r,,'^][r^/]) (31) 

Notice that, for the sake of simplicity, we have introduced a slight abuse 
of notation, since in [1] and [16] the same symbol -R^ instead denotes the 
inner curvature defined with the help of projections. Actually anything goes 
like in [1] section 4 with the Gauss-Codazzi framework, with the difference 
that here we don't have to make projections, which would involve product 
times a discontinuous tangent metric. Moreover here we don't even have 
to distinguish between the cases of E timelike or lightlike. In other words 
our mean-value differential geometry on a hypersurface is a very simple, in 
conceptual terms, analogue of the Gauss-Codazzi apparatus. 

Thus, with the heuristic theory of [1] section 6 (see also [1] for the timelike 
case) in mind as a prototype, we expect the jump of the Christoffel symbols to 
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play the main role, in place of the secund fundamental form, in the definition 
of compatibility conditions for very weak solutions of the Einstein equations. 
Indeed, this happens, as it will be shown in the following. 

4 Complex mean-value formalism 

The metric being dicontinuous on S, we are missing the fundamental tool 
to rise and lower indices, and to construct curvature in the traditional way. 
This is the reason why sometimes one is tempted to introduce some hybrid 
metric object on S to replace the metric, even in the (C°, piecewise C^) case 
(see e.g. |5j). It is reassuring to find out that the framework of the preceeding 
section can be confirmed by such a kind of approach. 

It would be desiderable to simply replace g with ^ on E, but it is easy to 
check that 'g has not the necessary algebraic requisites; in particular we have 
9af39°''' ^p'' ■ Consider instead: 

^a/3=^„/3 + ^(l/2)W, = r^ -i{l/2)[g^P] (32) 

where i is the imaginary unit (i.e. we have = —1). It is easy to check, 
with the help of ([3]), that we have: 

~9apr' = s^'+t[9ap\r' (33) 

i.e., in particular: '^{gapg"''') = Sf^^. For the sake of brevity in the following 
we will denote A ^ B the relation 3?(A) = 3^(5). Thus the pair ga/3 and 
g"^ is a good candidate replacement for the metric on E, for the purposes of 
rising and lowering indices. Now, similar to (132|) let us introduce: 

Tafsu = T^p, + z(i/2)[r„^,], f = r„^- - ^(i/2)[r„^-] (34) 

so that we have: ~ ^a^ug"^^ and conversely: Tapv ~ ^ap'^gua- Let us 

now introduce the differential operator V on S, which makes use of F in 
place of F. As we could expect we have: 

Vp^a;3~0, Vp^'^-O (35) 

which is the replacement on S for the covariant conservation of the metric 
tensor. 
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Now let us construct on E the complex curvature tensor R in the familiar 
way, but with F in place of the ordinary Christoffel symbols (which are 
undefined on E): 

RaPp" = djsTap'' — daTi3p'^ + T/j^^Tcp'^ — Tap'^TiSp^ (36) 

We rather unespectedly find out that 

^ = (i?s)./3p^ + ^(l/2)[i?„;3p'^] (37) 

i.e. in particular we have: Rapp" ~ {RT.)ai3p'^, where Rs is given by (131]) . 
This is just another reason for identifying i?s as the replacement for the 
curvature tensor of S, which is the first step of our path to the generalized 
compatibility conditions. 

5 Generalized compatibility conditions 

Let us now consider limit / ^ 0+ of the curvature tensor of the subdomain 
by ([2]) we have: 

{RaPp"")^ = Ra/Sp"^ + (1/2)[-Rq/3p'^] (38) 

and, by ([27D: 

{RaPp^y = RafSp" + V[;3[r,]/] (39) 

We also have, by ([SI]): 

(/2„^/)+ = (Rj:U,'^ + V[^[r,]/] + [r,[/] [r,]/] (4o) 

We see that R and i?^ only differ by terms proportional to [F], and not 
involving derivatives of it. Thus, in view of neglecting these tems, in the 
following we will consider R instead of R-^; this simply avoids the introduction 
of the symbol " = " ^ with the meaning of equality but for terms not involving 
derivatives of [F] (which here replaces the second fundamental form K) as in 
P section 6. Then for the Ricci tensor R^p = Rapp" we have: 

(i?^,)+ = R^p + (1/2)V^(5/[F,/] - [F^/]) (41) 
and for the curvature scalar R = Ra°': 
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R+ = R+ (l/2)V^([r/'^] - [r/^']) (42) 

Now, to construct the Einstein tensor we have to remember that, since 
the metric is also discontinuous: 

{9a^V = 9af3 + {l/2)M (43) 

so that we have: 

{G,,y = G,, + (1/2)V47^,/ - (1/8) W([r/1 - [r/^])} (44) 
where we have denoted, for the sake of brevity: 

n,/ = {5/[r./] - [r,/] - (i/2)^,^([r/i - [r/^]) } (45) 

Let us fix a coordinate chart and consider a generic (for the moment) regular 
prolongation for G, so that its mean value is defined in the whole Q. Now 
consider the Riemann 4- volume integral of over the domain from 
the Green theorem we obtain (for the general definition of integral on a 
hypersurface see |6] p. 6): 



X, G,, = G,, + (1/2) - (1/8) t^[g,,] {[Vr] - [r/^) 

(46) 

The analogous formula for Vt involves —I as the outgoing normal vector 
and the metric g~f^ = g^/j — (1/2) [(7^/3], so we have: 



l_ Gp, = l^^G,, + (1/2) Jj;n,,'^ + (1/8) JjM{[rn - [r/^]) 

(47) 

and consequently we have: 



„''* = 4.«^' + i.W (48) 

Thus reasons similar to those of the heuristic theory (see H section 

6) lead to the reasonable hypothesis that G remain bounded in the neigh- 
bourhood of S, for any admissible prolongation, so that from the volume 
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integral of the Einstein equations, with the presence of an eventual source 
term concentrated on E: 



/ Gaf3 — —X I Tap — X I Tai3 



(49) 



where x denotes the gravitational constant, we conclude that 



X 



(50) 



which is our heuristic reason for considering the following set of general- 
ized compatibility conditions to hold on S as a replacement for the Einstein 
equations: 



In the simpler case £q, G C°, it is very easy to check that the object if^HfSp'^ 
is gauge-invariant in the sense of (ITTl) . as we could hope. 

Turning now to the comparison with the C° case, we see from eq.s (71) 
and (85) of [Ij that our generalized conditions flSTl) are formally identical to 
ordinary compatibility conditions [eq. (110) of the same paper], if expressed 
in terms of [T] (which in the general case is a function of the jump of the 
metric [g] as well as of its weak discontinuity dg). Therefore it is clear 
that generalized compatibility conditions reduce to ordinary ones in case the 
metric is continuous, i.e. in case [ga/s] = 0. 

In particular, let us suppose g g(C°, piecewise C^) and / G C°(fi); let 
us moreover suppose {i ■ i) > 0, i.e. S timelike. By definition of Christoffel 
symbols, and from (11) of [Tj, we have: 



Q is a vector which can be set to zero with a suitable gauge choice; it plays 
no role in flSTl) . as one would expect, in fact we have: 




(51) 




gu'' + {i-i){Q-N) 
-gu'' + ii-i)iQ-N) 



(53) 
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and, since g = g = h{N) + N ® N, we have from fl45l) : 

if^Hf,/ = Gpp - h{N)ppGJ' (54) 
i.e., according to (88) of [1]: 

^tiHpp' = 'Hpp (55) 

as expected. Now let us instead suppose {(- ■ tj = 0, i.e. S lightlike. Let 
u G C^{Q) be a given auxiliary reference frame. From eq.s (21) and (16) of 
[Ij we have: 

and consequently, from (18) and (19) of the same paper: 

£^[r^/] = LpB{u,n)p + LpBiu,n)^-iu-ifLpLp{Q-L) 

£^[r,/] = Lf,Lpg{u,n),'^-{u-iYLpLp{Q-L) (57) 

We therefore have: 

= Gi^^ n),"LpLp - LpB{u, n)p - LpB{u, n)p (58) 

i.e. again, according to (83) of pLj, we have: dpTipp^ = Tipp, as expected. 

Therefore the set flSTl) of compatibility conditions, together with ordinary 
Einstein Equations to hold on each side of the discontinuity hypersurface, 
defines the class of generalized regularly discontinuous solutions of the Ein- 
stein equations. And in case [g] = 0, i.e. for continuous metric, from such 
conditions we recover the ordinary compatibility conditions for regularly dis- 
continuous weak solutions. 

However, in the generic case we have some differences, as we are going to 
show in the following. 

6 A class of spherical boundary layers 

Let us consider the match of two piecewise-C° regularly discontinuous spher- 
ical solutions of the vacuum Einstein equations, of the form 

ds'^ = -a(r, t)dt'^ + b{r, t)dr^ + r^dQ^ (59) 
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with dQ^ = dO"^ + sin^ Odip^, across a spherical admissible gravitational dis- 
continuity hypersurface S of equation r = p{t), with p{t) G C^. Therefore 
the form ^a = — P^a is continuous (while = g^'^la in general is not). 
We suppose globally C° coordinates, the same form of the metric in both 
domains fi"*" and and the identification = t~, r+ = r~, 6^ = 6~ , 
ip^ = on S. Leta, 6 > and let a,b E piecewise-C° be regularly dis- 
continuous on S and with regularly discontinuous first derivatives. Let us 
denote by a dot the partial derivative with respect to t, and by a prime that 
with respect to r. Let moreover condition a — bp > 0, i.e. {i ■ i) > 0, hold on 
both sides on S. 
We have: 

[9ai3] = -[a%%'+[b]5J5/,' (60) 

Now let us define the match as a generalized regularly discontinuous solution 
by (1511) . with T = 0, i.e. in the absence of stress-energy concentrated on E. 
In this case our compatibility conditions reduce to: 

i,[r,/]-QT,,^] = (61) 

which, for a match of metrics of the kind fl59|) . are equivalent to the following 
system: 

p[bb-^] + [a'b-^] = 
p[ba-^] + [a'a-i] = 
p[a'a~^] + [da~^] = 
p[b'b-^] + [bb-^] = 

i.e. we have [b] = and consequently: 

p[b] + [a'] = 
p[ba-^] + [a'a-i] = 
p[a'a^^] + [da""^] = 

p[b'] + {b] = 

and from ([3]): 

_ m + K] 

_ _ {pb + a')[a^ 

{p 17 + d)[a-^] + {p[a'] + [d])a-i 

m + [b] 



(62) 



(63) 











(64) 
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Now if we had both p[6] + [a'] = and p6 + a' = 0, by ([2]) we would have 
p6 + a' = on both sides of the hypersurface. We discard for the moment 
this singular situation, and from ( 16^ 9 we conclude that [a~^] = 0. 

Thus in this case our generalized compatibility conditions imply [a] = 
[b] = 0, i.e. they force the match to be piecewise-C^. 

In [Ij we have already studied some examples of piecewise-C^ matchs 
of metrics of the kind (j59l) at a hypersurface of constant radius r = ri,, with 
•^a = Sa^. Namely, we have considered: external Schwarzschild - internal 
Schwarzschild; external Schwarzschild - Tolman VI; external Schwarzschild - 
Tolman V. Such matchs obviously have ia ^ C^] moreover condition p b+a' ^ 
reduce in this case to a' 7^ 0, which is obviously satisfied. In each case we 
have verified that condition [a] = [6] = imply (9a = (where d denotes first 
order discontinuity), which then define the match as a boundary layer ^ (it 
actually also imply 96 = 0, as one can verify). Such is a general result, since 
for a metric of the kind (l59l) the completely temporal and radial components 
of the Einstein tensor are independent from the second derivatives of the 
metric: 



(65) 



Gtt = -a{b'r + b'^-b)/r^b^ 
Grr = —{a'r — ab + a)/ar 

so that the corresponding vacuum Einstein equations reduce to: 

b'r + b^-b = 0, 
a'r — ab + a = 0. 

Now, since in the match of (l59l) vacuum solutions equations (166!) are satisfied 
on each side of the interface E, their jump is in particular null, and from ([3]) 
we have: 



(67) 



n[b'] + (26 - 1)J6] = 
n[a'] -a[b] - (6+ l)[a] = 

from which it clearly follows that conditions [a] = [6] = imply [a'] = [b'] = 0, 
i.e. da = db = 0. 

Summarizing, for the match of two piecewise-C*^ regularly discontinuous 
spherical solutions, in the above hypothesis, generalized compatibility con- 
ditions (I5TI) imply [a] = [6] = i.e. they force the match to be C°. On 
the other hand conditions [a] = [6] = imply that E is a boundary layer. 
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Therefore for such spherical matchs generahzed compatibihty conditions flSlI) 
are necessary and sufficient for the match to be a boundary layer. 

7 Generalized gravitational shock waves 

Let us consider the match of two plane wave metrics of the form 

ds^ = -2d^d7] + F{ifdx^ + G{ifdy^ (68) 

across a hypersurface S of equation ^ = 0. Here ^ and rj are the two null 
coordinates. We suppose continuously matching coordinates and F, G reg- 
ularly discontinuous, together with their first and second derivatives. The 
gradient vector of S is the continuous characteristic (on each side of S) vector 
P = 5 ^ 

Generalized compatibility conditions flHT]) in the case T = (i.e. no stress- 
energy concentrated on the hypersurface) reduce to the following single scalar 
equation: 

[F-^F' + G-^G'] = (69) 

which characterize the generalized gravitational shock wave. Let us now 
study compatibility of (jHlD with the Einstein Equations. Einstein vacuum 
equations also reduce to a single scalar equation: 

F'^F" + G-^G" = (70) 

which we suppose to hold on each side of the hypersurface S; thus replacing 
F"*" and G^ by their expressions in terms of F, [F], G and [G] according to 
([2]) gives rise to the following couple scalar conditions: 

{2W+ [F"]){2G+[G]) + (2G^+ [G"]){2F + [F]) = (71) 
(2F^- [F"])(2G- [G]) + (2G^- [G'"])(2F - [F]) = (72) 

Equations (!7Tl) -( 172l) are compatible with (J69l) . i.e. the three equations set 
can be solved algebraically with respect to F, [G] and to any member of the 
pair [F, G), and the solution is not necessarily trivial. 

Finally let us notice that, if the additional condition [F] = [G] = holds, 
i.e. if the solution is C°, condition reduces to F-^[F'] + G''^[G'] = i.e.: 
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F-^dF + G-^dG = (73) 

which is the analogous condition for the ordinary shock wave, according to 
[1] section 10.5. 

8 Slow generalized gravitational waves 

Let us start trying to match two vacuum solutions of the kind (!68|) across 
the timelike (on both sides) hypersurface S of equation ^ = C- Again we 
suppose continuously matching coordinates, F, G regularly discontinuous to- 
gether with their first and second derivatives, and T = 0. This times gener- 
alized compatibility conditions include (IMj) and the following two additional 
scalar conditions: 

[FF'] = 0, [GG'] =0 (74) 
i.e., in terms of F, [F], G and [G], according to ([3]): 



F[F'] + [F]F' = (75) 
G[G'] + [G]G^ = (76) 

It is easy to check that the system fl75l) - fl76l) is not compatible with fl7T|) - fl72|) . 
in the sense that the whole system does not admit non-trivial solutions for 
F, [F], G and [G]. 

On the other hand we have proved in section 6 that a wide class of gen- 
eralized spherical matchs at a hypersurface of constant radius necessarily 
degenerate to a G^ match. 

Other examples of degeneracy have not been included in the paper for 
the sake of brevity, but at least it seems to be a hard task to construct a non- 
trivial generalized match across a timelike (on each side) hypersurface, with 
no stress-energy content. Such difficulty is certainly not a proof that this is 
an impossible task, but it makes us wonder whether such a solution should 
necessarily degenerate to a boundary layer, just like it happens for ordinary 
C° solutions (see e.g. [Ij). This would forbid the existence of generahzed 
solutions which propagate at a speed slower than light. Such would be a 
desiderable prohibition under certain respect, since one could expect that 
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gravitational interactions in vacuo must necessarily propagate at the speed 
of light also in a generalized theory. 

In general terms, since for generalized solutions the metric is discontinu- 
ous, a hypersurface can in principle have different signatures on the different 
sides; for this reason we cannot simply distinguish between the timelike and 
the light like for usual C° solutions. We should rather distinguish 

between three cases: timelike-timelike, timelike-lightlike (or conversely) and 
lightlike-lightlike. 

In any case it is legitimate to expect that, at least in the timelike-timelike 
case, similar to the timelike case of (C°, piecewise C^) solutions, absence of 
stress-energy concentrated on S should imply the solution to degenerate to 
a boundary layer p. 

Unfortunately for generalized solutions we still have no proof that absence 
of stress energy concentrated on S does necessarily imply the degeneracy of 
the solution to a boundary layer. 

Therefore, although the examples considered in this paper seem to sug- 
gest that such property could hold true also in the generalized case, for the 
moment such result is still a conjecture; we thus have to admit that the the- 
ory in principle allows propagation of generalized gravitational shock waves 
at lower speed than the speed of light. We would call such waves "slow gen- 
eralized gravitational shock waves". It would be reasonable to forbid this 
situation as unphysical, but for now this can only be done ad hoc, by means 
of a corresponding additional hypothesis. 

9 Non-symmetric stress-energy 

Notice that J^Tipp^ is not necessarily symmetric; from identity: 

= {l/2)g-'d^g (77) 
where g denotes the determinant of the contravariant metric, we have: 

I.nyppf = {l/A){lp[g-'dpg] - Ip[g-'dpg]) (78) 

Thus the generalized scheme allows in principle the presence of non sym- 
metric stress-energy on the discontinuity hypersurface. We will display non- 
trivial examples of non-symmetry in the following section. Notice that the 
right hand side of fl78|) is identically null in case g E and ia ^ C^, since 
in this case we have [g~^d^g] = i^g'^dg. 
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A non-symmetric Einstein tensor is a feature of Einstein-Cartan theory of 
gravitation (see [IH], see also [3j section 7.2), where it is due to the presence 
of torsion in the non-symmetric connection used to construct generahzed 
curvature. Thus the generahzed theory can be interpreted, at least to some 
extent, as introducing a torsion equivalent tool on the shell only, even if there 
are no geometrical objects in our theory which can be directly interpretated 
as torsion. However, Einstein-Cartan theory also has a spin - angular mo- 
mentum field equation in addition to the Einstein equations, which here is 
missing. 

In the literature, compatibility conditions for C° solutions of boundary 
layers [20], and recently of shock waves and thin shells [5T], have been studied 
also in the framework of Einstein-Cartan theory; actually this can lead to 
non-symmetric stress-energy on the shell. But in that theory this feature 
is inherited from the ambient spacetime, which is not here: non-symmetric 
stress-energy arises on the shell only, in consequence of the theory. This 
interesting feauture probably is worth investigating. 

10 Generalized thin shells 

Now let us consider a more general form of the spherical metric: 

ds"^ = -a{r, t)dt^ + 6(r, t)dr^ + c(r, t)dn^ (79) 

Let us consider a match of two spherical solutions of the Einstein equations 
across a timelike (on each side) hypersurface of equation r = p[t). Again we 
suppose pit) G and therefore £q, = 5^ — p6j G C°. Let the coordinates 
be globally, and let the metric have the same form (!79|) in both domains 
Q'^ and , with the identification = t~ , = r~ , 6^ = 6^ , ip^ = (f~ on 
S. 

Let moreover a,b,c > and let a,b,c G piecewise-C° be regularly dis- 
continuous on E and with regularly discontinuous first derivatives. Again we 
denote by a dot the partial derivative with respect to t, and by a prime that 
with respect to r. 

In this case for the left hand side of the generalized compatibility condi- 
tions i^TC^p^ we obtain: 
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+ (Ka-V2 + c'c-i] + p[6a-V2]) VV+ 

+ ([da-V2 + cc-i] + p[a'a-V2])VV+ ^on^ 

-([66-72] + p[6'6-V2 + c'c-i])VV+ 

+ ([c'6-V2] + p[ca-V2])(5/5/ + sin^ + 

-{\h-\a'a-^/2 + c'c-i)] + p[a-^(66-V2 + cc"^)])^^^ 
where, obviously: 

= -a5/(5/ + M/^p-^ + c{5p%; + sin^ ^5/5/) (81) 

We now mean to interpret flHOl) as the matter-energy of a thin shell. Let 
us first get back to the particular case p = (static shell) and a = h = 
c = 0, to make the interpretation simpler by eliminating the non-symmetric 
component; rearranging some terms we in fact obtain: 

^^.n^p^ = (-K&-V2] + ac-Hc'&:V2])W+ 

+ ([a'a-V2 + c'c-i] - hc-^[c'h-^ /2])5p'5;+ (82) 
+ (c-i[c'6~V2] - [h-\a'a-'/2 + c'c"!)])^^. 

This can be interpretated as a perfect isotropic magneto-fiuid thin shell with 
infinite conductivity, i.e. we can solve the compatibility conditions by con- 
sidering the following stress-energy as the right hand side: 

Tap = {po+P + ph^)UJJp + {l/2)^h^)gap - pKhp (83) 

where po is the proper density, h the magnetic field and p the magnetic 
permeability [221 [23l EH ESI E]; here we define = Kh^g'^'^ . In fact it 
suffices to define the following 4-velocity vector: 

Ua = V^a- (l/4)[a2]a~i5„* = (^f/^dj (84) 

which by construction is unitary on S, in the following sense: UaU/s'g'^^ = — 1, 
and the following magneto-hydrodynamical variables: 
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= x'^^[a'b-^]/2_+ X'^c-^b FT/2 - aFT + l)[c'b~^]/2+ 
-X~M^~^](Q^/2 + c^) - (3/2)x"^FT[a'a-V2 + c'c-^] 
p = x'^c-\bF\/2 - l)[c'lr\/2] + (l/2)x"^FT[a'a-V2 + c'c-^]+ (85) 

;i„= ±v^FTx-HKa"V2 + c'c-i] - 5c-i[c'6-V2])5„^ 

to match ^ and ([Ml) via f^H/jp'^ = -x^^p- If [«] = [b] = [c] = then the 
generahzed shell fl85|) degenerates to the C° magnetohydrodynamical shell 
considered in jl] section 10.1, in the particular case p = 0. 

The slightly more general case ofd = 6 = c = 0, but p 7^ 0, displays 
non-symmetric terms in flHOj) : however it is not difficult to see that the per- 
fect magnetofluid interpretation still holds, provided such additional non- 
symmetric terms are interpreted, or neglected. In fact in this case we have: 



\,n^/ = (-K6-V2] + a c-i[c'rV2]) W+ 

+ {[a'a'^/2 + c'c-i] - 6c-i[c'rV2])5^^'5/+ 
+ (l/2)p[a'a"V2 - b'b-'/2 - c'c-i](5^^'V + 
+ (l/2)p[a'a-V2 + b'b-'/2 + c'c-i](5/'V - 
+ (c-Hc'6-V2] - [rHa'a-V2 + c'c'')])g^^ 



) + 
) + 



(86) 



Now let us consider, for the sake of brevity, the following quantities: 



a 



P 



4" ha 2b 



mi] - [|])^«-^ 



-\-] 

c I2bi 



hbi 



2a c 



c' 
26 



1 \— 
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2b 



c'l2 



(87) 



(88) 



cL26J L26J 

and let us suppose that inequality a < holds, which is necessary for the 
physical interpretation. In fact in this case the following vector: 



■ c L2feJ 



-p\- 
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2b 



is a unit timelike vector on S, in the sense that UaUpg'^^ 



(89) 

-1. Rearranging 



terms, fl86D now reads: 
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l.-Hp,^ = aU,U, + P6,%^ + \p[i + ^] (5/5/ - 5,*5/)+ ^^^^ 

)9l3p 



which can be matched via i^Hisp^ = —xTpp with a stress-energy tensor of 
the following kind: 



Tap = {po+P + fxh'^)UaU(3 + {p + {I / 2) ^h'^) Q - fxhahp + (91) 
where A denotes the anti-symmetric term. We have: 



PO - A " I A -A2bi A I" \2a ' c JJ_ 2x 



P = -X-'m + X-'[b-\i + i)] - ^f3b~' (92) 



while the anti-symmetric term A reads: 



1 



a' b' c' 



The interpretation of such term is still missing; alternatively it could be 
neglegted by adding the further hypothesys: 



a' h' c'- 
2^ ^ 26 ^ 7. 



(94) 
which is equivalent to [g^^g'] = 0, as we could expect from (178!) . 
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